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Abstract - Oppositely driven binary particles with repulsive interactions on the square lattice 
are investigated at the zero-temperature limit. Two classes of steady states related to stuck config- 
urations and lane formations have been constructed in systematic ways under certain conditions. 
A mean-field type analysis carried out using a percolation problem based on the constructed 
steady states provides an estimation of the phase diagram, which is qualitatively consistent with 
numerical simulations. Further, finite size effects in terms of lane formations are discussed. 



Introduction. — Noncquilibrium systems exhibit col- 
lective phenomena under certain conditions. The under- 
standing of the relationship between such collective phe- 
nomena and microscopic properties of the system has been 
an important topic in statistical physics. Recently, it has 
been observed that various systems consisting of oppo- 
sitely driven binary particles with repulsive interactions 
exhibit a collective phenomenon called lane formation, 
where the same type of particles align to the driven di- 
rection [H7] . The observation of lane formations in vari- 
ous systems suggests that certain types of properties are 
universal in lane formations, irrespective of the systems. 

Thus far, a phenomcnological theory for determining 
a special condition to show lane formations in colloidal 
suspensions has been presented, which is consistent to nu- 
merical results in some parameter regions (2JI31 ■ However, 
we still lack the knowledge for the relationship between 
lane formations and microscopic properties. For exam- 
ples, questions on whether the special condition is related 
to some nonequilibrium phase transitions or which types 
of fluctuations, for instance, finite size effects, appear in 
lane formations remain unanswered. 

As is common in the literature on statistical physics, 
lattice models are often useful for obtaining insights into 
universal phenomena, owing to the simplicity of the sys- 
tem and the universality itself. For example, driven lattice 
gas models and simple exclusion processes have provided 
many insights into universal aspects of nonequilibrium sys- 
tems, such as phase transitions and the KPZ universality 
TfJl . In this context, a simple lattice model exhibiting a 



lane formation can be useful for understanding lane forma- 
tions more comprehensively. However, extensive studies 
have not been conducted on this topic. 

In this paper, in order to obtain insights into the univer- 
sal aspects of lane formations, we propose a lattice model 
where binary particles with purely repulsive interactions 
are driven in opposite directions. We present the exact 
constructions of two classes of steady states related to 
stuck configurations and lane formations. On the basis of 
these constructions, we introduce a percolation problem in 
a stochastic cellular automaton, which helps us estimate 
the phase diagram. Further, we discuss finite size effects 
in lane formation using the constructed steady states. 

Model. — Let us consider a square lattice A consisting 
of site i G {{i x , iy} G N x N|l < i x , i y < L} where N is the 
set of natural numbers. In order to introduce binary (pos- 
itive and negative) particles on the lattice, we consider an 
occupation variable <7j G Z (— L 2 < er.; < L 2 ), where Z is 
the set of integers. In this expression, Ui indicates that 
there are |<7j| number of particles with the sign of ai at 
site i, and no particles if <Ji = as illustrated in figure [TJ 
Therefore, the density p of the particles in the system is de- 
fined as p = J2ieA \ a i\- We consider a situation where 
each particle interacts with soft-core repulsive interactions 
with the following Hamiltonian for a = {(Ji}i<£A'- 



ff(ff)=v y>i-i(k<i-i), 



(1) 
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Fig. 1: Schematic of the model showing the correspondence 
between the configuration of {(Ji}i (left) and that of binary 
particles (right). 

where Vo > 0. On the basis of this purely repulsive 
Hamiltonian, we consider continuous-time Markov pro- 
cesses where a particle hops to another site with some 
transition rates; this process is expressed as follows. Pre- 
liminarily, in order to express the hopping of one particle 
from site i to site j, we define an operator Fij that sat- 
isfies FijGi = Gi - sgn(erj), FijCFj = Gj + sgn(<7i), and 
FijfJk = where sgn(cc) = — 1 for x < and sgn(x') = 1 
for x > 0, otherwise sgn(x) = 0. Here, the destination 
site j can be in set B^ = {j e A||i — j\ < /hop} where 
\i - j\ ee y/(i x - j x ) 2 + (i y - j y ) 2 and / hop > 1. as shown 
in figure [2l The transition rate R(cr — > F^er) between 
state er and Fijcr with cr, ^ is 

Q(aiaj)wij(ai,aj)r((T -t i*y<x), (2) 

with r(er — > Fijcr) = 

exp(-/3/2(Afly(«r) - D^) - AE^cr))), (3) 
Wij (at ,aj) = dij J^i±2d ) (4) 

where Ai^(er) ee (H(F tJ cr) - H{a))/\i - j\, Dyfa) = 
f(jx ~ *a;) s g n ( cr i)/|* ~ ill and 6 is a step function such 
that 0(x) = 1 for x > 0, otherwise 0(x) = 0. We define 
r(er — > Fijtr) = with at = to maintain consistency 
in this description. Further, we assume that AE^er) = 
AEj i (Fija) holds; an explicit form of this expression is 
provided later. According to these properties, one can 
prove that R(cr — > F^cr) satisfies the detailed balance 
condition with / = 0. 

Here, we explain the physical meaning of each term. 
First, the term Dij(<Ji) governs the tendency of positive 
(negative) to hop in the positive (negative) direction along 
the a;-axis with a driving field /. Second, the weight 
Wij{(Ji, (jj) along the edge for sites i and j has two phys- 
ically reasonable implications. One implication is to as- 
sume that the diffusion of a free particle is uniform for the 
x(y) and diagonal directions. Explicitly, we can realize 
this situation by considering that \i — j\dij is independent 
of j, for example, using X^eB K — il^j = ^- wmcn leads to 
= 1 with /hop = 1 or d^ = \i — 7 1 1 / 2 with /hop = \/2. 
It should be noted that dij = dji by definition. The other 
implication is that the term |<jj +<Jj \ ensures that the num- 
ber of hopping events at the edge for sites i and j occur 
in proportion with the number of particles along the edge. 
Third, the step function Q{uiUj) ensures that every type 
of particle is conserved. 
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Fig. 2: Hopping from i to j with Z hop = V2: AE ( 'j(a) is calcu- 
lated by AEij (a),AE i k 1 (er), AE %k2 (cr), AE ik3 (tr), AE iki (er), 
AEjiiFijtr), AE jkl (FijCT), AE jk2 {F ijC T), AE jk3 (F ijC r), and 
AE jk4 (F ij( r) (left), or AE t] (a), AE ikl ((r), AE ik2 (cr), 
AEjiiFijO-), AE jkl (F ijt r), and AE jk2 (F ijC r) (right). 

Here, we consider an explicit form of Ai?° (cr). In order 
to construct a model with universal properties, we develop 
a guiding principle for determining forms of AS?-(cr) in 
the following manner. The idea is to consider the systems 
described by overdamped Langcvin equations where parti- 
cles move in the exact direction of the driving force at the 
infinite driving force limit pQ. Keeping this in mind, we 
assume that particles at each site i hop only to sites fc 6 Bj 
with maximum Dik(<Ji)@(<7iCrk) at the infinite driving field 
limit, which we refer to as normality of the limit. For ex- 
ample, the following expression provides the normality of 
the limit; 

AE°j{a) = mm {AE ik (a),AE jk {FijCT)}, (5) 

where AE ik (cr) = AH ik (<r) - D ik ((Ti), and = {k E 
A\\k — i\, \k — j\ < /hop} with /h op = This ex- 

pression also satisfies A£° (er) = AE^i^er) < and 
< r(cr — > Fij<j) < 1. In this study, we investigate 
only this case. In order to capture a physical image of 
AEfj(a), we show AE ik required to calculate AEf (cr) in 
figure [2j It should be noted that in order to ensure that 
only AE°j(cr) = AE^F^cr) holds, we can set A£° (er) as 
a constant. However, in this case, the model does not fol- 
low the normality of the limit. As a result, the arguments 
mentioned in the rest of this paper do not hold. On the 
other hand, we expect that the arguments mentioned in 
the rest of this paper are robust against small changes in 
the expression for AEfj (a) if the expression provides the 
normality of the limit] however the extent of changes for 
which the argument hold is uncertain. This type of guiding 
principle used to determine transition rates with universal 
properties has been discussed in different nonequilibrium 
lattice models with a nonconscrved variable 

In this study, we consider the behaviours of the system 
by changing the set of parameters (p, /). Further, for sim- 
plicity, we focus on the case with the zero-temperature 
limit, the periodic boundary condition, and 50:50 binary 
mixtures where the densities of positive and negative par- 
ticles are p/2. If / = 0, the steady states are determined 
by the canonical distribution of Hamiltonian ([T]). In par- 
ticular, one can immediately determine that if / = and 
P < 1, o~i in equilibrium states randomly takes only the val- 
ues of 1, —1, or with a fixed density p. Obviously, there 
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Fig. 3: Particle configuration of a numerically observed global- 
MC steady state at / = 100.0 and p = 0.5 with L = 50. 
The red (green) circles are positive (negative) values for the 
occupation variable. 

are no singular points in the parameter space under this 
equilibrium condition. We consider relaxation behaviours 
from this equilibrium initial condition. It should be noted 
that these behaviours are still non-trivial even under these 
simple conditions owing to purely nonequilibrium effects, 
as explained later. While performing Monte Carlo simu- 
lations for this model, we randomly select a particle. If it 
is located at site i, we select any j G B^ with probability 
dij/(2 + a/2)- Then, both the directions from i to j and 
from j to i are adopted with a probability of 1/2. Finally, 
a particle hops to the adopted direction (if j — > i) with the 
probability 0(<7i<7j)r(<T — > FjiCr), and if there are no par- 
ticles at site j, nothing occurs. This process is repeated 
and time t ~ 1 corresponds to the repeated L 2 p/r a steps 
where r = 1/(4 + 2^2). 

Exact arguments and the related conjectures. 

Let us confirm some non-trivial arguments for the model. 
Preliminarily, we consider a current-like quantity as fol- 
lows. Let us define a colored current J[t) at time t with 
<r as 

J ( i ) = Tl _ 5Z ®(o-zO-j)sgn(o-i)r(<r F^ct), (6) 

where — 1 < J(t) < 1 by definition. We consider J = J(t) 
such that J(t) = lim To ^oo ^- Y^vZt ^(0- ^ i s trivial that 
if / = 0, J should be absolutely zero with r — > oo owing 
to the symmetry of the directions. By using this quan- 
tity, we define the maximum-current (MC) steady states 
as J = 1 and the stuck steady states as J — 0. Next, 
using the colored current, we construct the microscopic 
characterizations of such steady states in some parame- 
ter regions. After these constructions, we present some 
conjectures related to the constructed steady states. 

One explicit class of stuck steady states with < / < 
2Vb for p < 1. Let us construct one explicit class of the 
stuck steady states at < / < 2Vq in the following man- 
ner. We consider an auxiliary variable bi E { — 1,0,1,2} 
as follows. 

{fri = 2} = {o- ix4y = 1, o- ix+ i,i v = -1}, (7) 
{N < 1} = Wi = (8) 



Fig. 4: Particle configuration in a numerically observed stuck 
steady state at / = 1.0 and p — 0.5 with L = 50. 

On the basis of this variable {bi}i, we consider a cellular 
automaton for determining the value of bi under the initial 
condition <5> = 2 and bi = for all the other sites i; the 
rule to update the values for all the other sites is explained 
in the following sentences. First, we set o'io_ lji o +1 = 2 or 
o>, i„+i = 2 or <7jo +1 10+1 = 2. Next, if b lxA a +1 = 2, we 

IB » H x 1 y y 

set a ix -i t ia +2 = 2 or b lx , iy+2 = 2 or 0-^+1^0+2 = 2. This 
process is repeated. After setting <Ti = 2 for at least one 
site for every values of i y , we set &i = —1 for i l y < i x < i r , 
where i l y = max a , i x + 1 such that bi =2 for a given value 
of i y and i T = max s i + 1. Similarly, we set <f , = 1 for 
i l < i x < i y ; where i y = max^ i x such that bi = 2 for 
a given value of i y and i l = max y i y . Thus, we obtain a 
sequence {bi]i. Finally, if we set H(cr) = 0, this system 
reaches a steady state with J = at the long time limit 
because {bi}i are independent of t in this construction. 
This corresponds to the fact that there exist the stuck 
steady states characterized by {bi]i with H(cr) = for 
/ < 2Vq and p = J2i \ a i\- Further, a state constructed 
by using {bi}i and {b'^i with b[bi = and H(cr) = can 
also be a stuck steady state. We define the stuck steady 
states constructed in this way as CM-stuck steady states. 
However, there are no principles to determine whether CA- 
stuck steady states appear as steady states under the equi- 
librium initial condition, which is discussed later. 

One explicit class of maximum- current steady states 
with sufficiently large values of f for any density. 
Let us construct one explicit class of the maximum- 
current steady states with sufficiently large values of / 
for any density, as follows. First, let us define N p , y = 

J2i> and N n , v = J2i' x Wi'\@(-<ri')$i' x ,y For 

eacfi value of y, we set N PtV and N ntV such that N PtV N nty = 
under the condition J2 y N p , y = J2 y Nn, y = pL' 2 /2. 
Then, if we set / such that 

V2 — 1 v 

we can easily find that this system reaches a steady state 
J = 1 at the long time limit because N PjV and N UiV are 
independent of t. This corresponds to the fact that the set 
{Np^y, -?Vn,2/}o<y<L under the condition mentioned above 
determines a maximum-current steady state for / satis- 
fying © and p = j,, J2y=o( N P,y + N n,y)- We define the 
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maximum-current steady state constructed in this way as 
global- MC steady states. However, there are no principles 
to determine whether global- MC steady states appear as 
steady states under the equilibrium condition, which is 
also discussed later. 

Conjecture A: Mean-field type analysis of phase transi- 
tions. So far, we have microscopically constructed the 
CA-stuck steady states and the global- MC steady states, 
which lead us to expect that singular points of J exist 
in parameter space (p, f) where the steady state changes 
qualitatively. With this background, let us assume that on 
fixing p, there is a line f c (p) where J is singular in driving 
field /. On the basis of this assumption, we attempt to 
approximately estimate f c {p) as follows. 

In low-density regions, particles appear to be driven in- 
dependently. Therefore, a non-zero value of f c (p) would 
be observed only in sufficiently high-density regions, and 
there should be a threshold density p t h below which 
/ c = 0. We estimate pth as follows. First, it should 
be noted that it is plausible for a type of seed to ex- 
ist under the initial condition for realizing the CA-stuck 
steady states. One reasonable candidate for such a seed 
is the particles present at the sites indicated by the one- 
step run of the cellular automaton from a focused parti- 
cle, which we refer to as an automaton-pointed segment 
of the focused particle, as discussed in the construction 
of the CA-stuck steady states. We assume that if there 
exists at least one sequence of automaton-pointed seg- 
ments starting from a focused particle, which connects 
with a site with distance L from the focused particle un- 
der the initial condition, the system can reach a CA-stuck 
steady states after the dynamical particle-exchange pro- 
cesses. As mentioned in the construction of the CA-stuck 
steady states, this consideration is applied in the case of 
/ < 2, which leads to f c {pth) = 2Vb- On this assump- 
tion, we consider the probability P y that there exists at 
least one sequence of automaton-pointed segments start- 
ing from a focused particle at i y = y, which connects with 
a site with i y = L under the initial condition. This is 
a typical percolation problem in a stochastic cellular au- 
tomaton with a percolation point p p |12j . Reminding that 
the number of sites in such an automaton-pointed seg- 
ment is 4, if each site in an automaton-pointed segment 
assumed to be independent, we can obtain the recursive 
equation P y = p(l — (1 — P y+ i) 4 ). Thus, we can esti- 
mate a lower bound p p = 1/4 of the percolation point 
p p , which also gives an approximate value of p t h- On 
the basis of this consideration, let us estimate / c (p) at 
p > pth- Here, we consider the condition with which 
the percolation breaks owing to the hopping of a parti- 
cle through an automaton-pointed segment, where only 
the same type of particles exist. Let 5n denote the incre- 
ment in the number of particles in an automaton-pointed 
segment when p is increased by 5p from p t h- Then, we 
can estimate Sn = 48p approximately where 4 is the 
number of sites in one automaton-pointed segment. In- 



deed, if / is increased by Sf from / = / c (pth) satisfying 
Sf > V (n + 8n)(n + Sn - 1) - n(n - 1)) = V (nSn + 8n 2 ), 
where n = 1 is the maximum particle number per one site 
at / < /c(/3th): a particle can pass through an automaton- 
pointed segment in the driven direction even for the worst 
case. Thus, we obtain f c {p) ~ /mf(/°) where 

/mf(p) = W (p - p\) + 16Vb(p - p\ f + /c(pth), (10) 

for p > p l p , otherwise /mf(p) = 0. 

Further, using ([9]), we can obtain the minimum possi- 
ble value of / = f = ~ 6.83 to realize MC-global 
steady states at the dilute limit. At / > f . microscopic 
events change considerably, where a particle can overlap 
with the front particle in the driven direction without be- 
ing affected by other particles. As a result, fluctuations 
in the local density could increased considerably, possibly 
causing another singularity in J even in finite density re- 
gions. Thus, we can predict that there is a cross point 
po ~ 0.69, determined by /mf(po) = jo, where the perco- 
lation effects and the large fluctuations of local densities 
become comparable. 

Conjecture B: Finite-size effects due to non- 
commutativity between t — > oo and L — > oo. Let 
us consider the situation with a sufficiently large driving 
field, where N p . y N niV ^ for some values of y at a finite 
time t, and focus on such a line-y, which is a set of sites 
{ixi 2/}i x • We consider the conditions where such all line-y 
can be divided into finite number of segments, each of 
which satisfies — where N®, s is the number 

of positive (negative) particles in segment x at i y = y. 
Clearly, there are no positive contributions to J(t) at 
the contact points between such segments. However, 
owing to the finite numbers of segments, we can neglect 
the effects brought about by such contact points with 
L — > oo. Hence, in the thermodynamic limit L — ¥ oo, 
J = 1 can be realized in other states besides global- MC 
steady states if we assume that the number of segments 
divided by L goes to zero at any time. Furthermore, 
when we consider the limit L — > oo at a fixed large value 
of t, such states appear not to relax to any global-MC 
steady states. This is because, in simple terms, the 
time required for the displacement of particles to form 
global-MC steady states would be at least of the order 
of L. Therefore, owing to the uniform configurations 
under the initial conditions, this type of transient state 
should be dominant for sufficiently large driving fields, 
as compared to the global-MC steady states, in the 
procedure where t — ¥ oo after L — > oo. 

Numerical simulations. — Next, in order to verify 
the plausibility of the obtained results, we perform numer- 
ical simulations with t = 50000to Monte Carlo steps and 
practically set r = 40000r . Further, we set Vq = 1 with- 
out loss of generality. We have verified that the following 
results are qualitatively unchanged with r = 20000ro. As 
shown in figure G2 we have observed a global-MC steady 
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Fig. 5: Colored current J averaged in the time period from 
t to 1.25r with L = 200. There is a strong sample-to-sample 
fluctuation from 2 < / < 4 with p = 0.6. 



state at p = 0.5, / = 100.0, and L — 50. It should be noted 
that observing global- MC steady states becomes consider- 
ably difficult for larger system sizes, which is discussed 
later. On the other hand, as shown in figure |H although 
we have observed a stuck steady state at p = 0.5, / = 1.0, 
and L = 50, the stuck steady state is not identical to 
the CVl-stuck steady states. Generally, some parts of CA- 
stuck steady states and other complicated configurations 
exist in the observed stuck steady state. Nevertheless, 
maximum-current steady states and stuck steady states 
exist in different parameter regions, which indicates the 
existence of phase transitions. 

In order to investigate the expected phase transitions 
explicitly, we change the external field by an increment 
of Sf = 0.5 at a fixed density p and measure AJ'(f) = 
J{f) — J{f ~ Sf). Then, we attempt to detect a point 



f n (p) = arg max, AJ(f) 



(11) 



which can be a reasonable candidate for the singular points 
f c (p) of J. It should be noted that the following observa- 
tion is made by one simulation run under an initial con- 
dition with L = 200. As shown in figure [SJ the numerical 
simulations demonstrate that it is easy to determine f n {p) 
at sufficiently low densities because J jumps discontinu- 
ously at some values of /. As shown in the figure |6l using 
f n (p), we can approximately estimate f c {p), which oper- 
ationally defines a flowing phase and a blocked phase. It 
should be noted that another weak jump around fo was 
observed independent of p at sufficiently low densities, as 
expected by the Conjecture A. However, it is not easy 
to estimate the location of f n (p) for p > 0.5 because of 
sample-to-sample fluctuations in f n (p)- For this region, 
we need another quantity for the plausible estimation of 
f n {p), as discussed later. Further, discontinuous jumps 
of J seem to disappear in high density regions such as 
p = 0.9 and 1.0. 

Next, we investigate finite size effects in the system, in 
particular, focusing on lane formations. Lane formations 
indicate that particles tend to be aligned to the direction of 
the driving field. The lane formations can be characterized 
by at least two approaches. The first approach, referred to 
as local lane formation, involves estimating how often the 



Fig. 6: Phase diagram estimated by / n (p) numerically (points) 
and by /mf(p) using the mean- field type analysis (real line). 
Observed weak singularities around fo in J and f\ are plotted 
(See the text for the definitions). The interval for fixed values 
of p is 0.1. Further, the data with p — 0.29 and p = 0.05 are 
plotted because/ n for p < 0.28 and the singularity around fo 
for p < 0.04 are not observed. 



same type of particles are nearest neighbors of the focused 
particle in the driven direction [7j . The second approach, 
referred to as global lane formation, involves estimating 
the number of times the same type of particles occur at 
all sites in the driven direction. [IH3] . As one possible way 
to quantify such local lane formations, we focus on 



$! = X (U - <h{lo), 



Oi' i ) and Y. 
»')• Here, 



where 0x('o) = 

Si' x ,i x )Q{±a ix ^ 
Si',i v )®{±<ri*,iv a 
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xf+x~ 



(12) 

and Y (M = 



with Xf 
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in order to determine 
nearest-neighbor particles for the focused particle, we set 
lo = 1 J for L/2 > [yo 1 J under the assumption that 
there two particles exist in the segment at an average dis- 
tance of approximately 21q in the driven direction. On the 
other hand, as one possible way to quantify global lane 
formations, we focus on 



% = ML/2). 



(13) 



As shown in figure [7l $i shows a clear discontinuous jump 
between a negative value and a positive value around the 
obtained f n (p) even at p > 0.5. Further, another singu- 
larity around fo is clearer than that of J . In high-density 
regions such as p = 0.9 or 1.0, it is also clear that there 
is one large discontinuous jump between a negative value 
and a positive value whose critical fields f\ are plotted in 
the figure |5] as a reference. More importantly, the global 
tendencies of $i do not depend on the system sizes, which 
indicates that behaviors of local lane formations are inde- 
pendent of the size of the system. We have numerically 
observed that global tendency of J also do not depend on 
the system sizes at / > f\(p), which is reasonable because 
J is also a local quantity. Therefore, we used f\{p) as 
the value of f n {p) for 0.5 < p < 0.8 in the thermodynamic 
limit. On the other hand, as shown in figureEl $ g shows a 
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Fig. 7: The local lane order parameter $i averaged in the time Fig. 8: The global lane order parameter $ g averaged in the time 
period from r to 1.25r for different system sizes at p = 0.3. period from r to 1.25r for different system sizes at p = 0.3. 



strong dependence on the system sizes. In particular, the 
tendency of these finite-size effects indicates that global 
lane formations disappear in the procedure where t — > oo 
after L — > oo, which is consistent with Conjecture B. 

As shown in figure |BJ /mf (p) and po obtained by the 
mean-field type analysis in Conjecture A provide qual- 
itatively consistent behaviours with f n (p), although the 
mean-field type analysis appears to overestimate the pa- 
rameter regions of the blocked phase. These deviations are 
rather natural for this mean-field type analysis because of 
the existence of spatial correlations that we have ignored. 

Concluding remarks. — We have proposed a simple 
lattice model for oppositely driven binary particles with 
purely repulsive interactions. We have exactly constructed 
global- MC steady states at sufficiently large values of the 
driving field for any density and CM-stuck steady states 
at small driving fields for low densities. A mean-field type 
analysis leads us to estimate singular points in the col- 
ored current, which are qualitatively consistent with nu- 
merical simulations. This strongly suggests that such a 
singularity originates from percolation in a stochastic cel- 
lular automaton buried in the equilibrium configuration. 
Further, we have presented a conjecture for the absence 
of the relaxation toward global- MC steady states (global 
lane formations) from the equilibrium initial conditions in 
the procedure where t — > oo after L — > oo, which is also 
consistent with numerical simulations. 

Here, we present some remarks on the universality of 
the observed phenomena in the proposed model. First, we 
describe the singular point /o in low-density regions, above 
which the driving field is large enough to allow a particle to 
enter inside a soft-core repulsive potential via a two-body 
collision. After such a strong collision, the fluctuations 
in the local density could be considerably strong. This 
type of strong density fluctuations might be related to in- 
stabilities of homogeneous density profiles in a Langevin 
system [3] and a two-lane traffic model [13] . On the other 
hand, the Langevin system exhibits a lane formation also 
in high-density regions [3], which might be related to the 
singularities above po in the proposed model. Further, 
in the case of adding a finite-range repulsive interaction 
S<ij> l^i cr i cr il to the Hamiltonian, we have numerically 



found that the lane formation is qualitatively similar to 
those in this study if \Vi\ <C f. Such extensions includ- 
ing finite temperatures are interesting when comparing the 
proposed model to the other more realistic systems with 
various effects such as inertial effects [33], attractive in- 
teractions [TS], hydrodynamic interactions [IB]. Thus, the 
robustness of the universality of the observed singularities 
in the proposed model are interesting topics for future 
studies. 
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